Numerical indications on the semiclassical limit of the flipped vertex 
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We introduce a technique for testing the semiclassical limit of a quantum gravity vertex amplitude. 
The technique is based on the propagation of a semiclassical wave packet. We apply this technique 
to the newly introduced "flipped" vertex in loop quantum gravity, in order to test the intertwiner 
dependence of the vertex. Under some drastic simplifications, we find very preliminary, but 
surprisingly good numerical evidence for the correct classical limit. 



Suppose you are explicitly given the propagation ker- 
nel Wt{x,y) of a one-dimensional nonrelativistic quan- 
tum system defined by a hamiltonian operator H 

Wt{x,y)^{x\e~^^'''\y) (1) 

and you want to study whether the classical {h — > 0) limit 
of this quantum theory yields a certain given classical 
evolution. One of the (many) ways of doing so is to 
propagate a wave packet ^px.pix) with Wt{x, y). Suppose 
that in the time interval t the classical theory evolves the 
initial position and momentum Xi,pi to the final values 
Xf,pf. Then you can consider a semiclassical wave packet 
4'xi,pi{y) centered on the initial values Xi,pi, compute its 
evolution under the kernel 

<P{x) J dy Wt{x,y) ^P.,,pM (2) 

and ask whether or not this state is a semiclassical wave 
packet centered on the correct final values Xf,pf. In this 
letter, we consider the possibility of using this method 
for exploring the semiclassical limit of the dynamics of 
nonperturbative quantum gravity. 

An explicit expression for the vertex amplitude in loop 
quantum gravity has been recently proposed in [T]; we 
call it here the "flipped" vertex, following [2]. This ex- 
pression corrects certain difficulties that have emerged 
with the Barrett-Crane vertex amplitude . In particu- 
lar, the Barrett-Crane vertex was shown to have a good 
semiclassical behavior as far as its dependence on spin 
variables was concerned, but not in its intertwiner sector. 
The flipped vertex W{jnm,in) is a function of ten spin 
variables jnm where n,m = 1,...,5 and five intertwiner 
variables i„ and it is hoped to correct the intertwiner 
dependence of the Barrett-Crane vertex. Here we inves- 
tigate this intertwiner dependence. A number of variants 
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of the flipped vertex have appeared in literature j4] , but 
we do not consider these variants here. 

The derivation of the vertex amplitude presented in 
[2\ indicates that the process described by one vertex 
can be seen as the dynamics of a single cell in a Regge 
triangulation of general relativity. This is a fortunate sit- 
uation, because it allows us to give a simple and direct 
geometrical interpretation to the dynamical variables en- 
tering the vertex amplitude, and a simple formulation of 
the dynamical equations. The boundary of a Regge cell 
is formed by five tetrahedra joined along all their faces, 
thus forming a closed space with the topology of a 3- 
sphere. Denote the area of the triangle (nm) that 

separates the tetrahedra n and m. Denote a^^''^^'' the 
angle between the triangles (mp) and {qr) in the tetrahe- 
dron n. And denote Qnm the angle between the normals 
to the tetrahedra n and m. These quantities determine 

entirely the intrinsic (^nm, ai™^'*'^'') and extrinsic (Qnm) 
classical geometry of the boundary surface. 

The ten spins jnm £^re the quantum numbers of the 
areas Anm [21 IS]- The five intertwiners i„ are the quan- 
tum numbers associated to the angles ai™^'''^^ More 
precisely, they are the eigenvalues of the quantity 

^{mp,qr) _ JS^^^ _|_ 2A„ipAqr COS Q!^'"^'''"-' . (3) 

Each tetrahedron has six such angles, of which only two 
are independent (at given values of the areas); but the 
two corresponding quantum operators do not commute 
[6J and a basis of the Hilbert space on which they act can 
be obtained by diagonalizing just a single arbitrary one 
among these angles. Therefore the intrinsic geometry 
of the boundary of a classical Regge cell is determined 
by twenty numbers, but the the corresponding quantum 
numbers are only fifteen: the fifteen quantities jnmi^n- 
These are the fifteen arguments of the vertex. When 
using the intertwiners i„, we have of course to specify to 
which pairing i^™^'*'^^ we are referring. 

The equations of motion of any dynamical system can 
be expressed as constraints on the set formed by the ini- 
tial, final and (if it is the case) boundary variables. For 
instance, in the case of the evolution of a free particle 



2 



in the time interval t, the equations of motion can be 
expressed as constraints on the set (xi,pi,Xf,pf). These 
constraints are of course m(xf — Xi)/t — pi = pf. (For 
the general logic of this approach to dynamics, see [8].) 
In general relativity, the Einstein equations can be seen 
as constraints on boundary variables Anm,cth^^''^^^ and 
Qnrn- Thcsc, iu fact, can be viewed as the ensemble of the 
initial, boundary and final data for a process happening 
inside the boundary 3-sphere. Such constraints are a bit 
difficult to write explicitly, but one solution is easy: the 
one that corresponds to flat space and to the boundary 
of a regular 4-simplex. This is given by all equal areas 
Anm — io, all equal angles i„ = zq, and 0„m = O, where 
elementary geometry gives 



V3 



Jo, 



cos O = 



(4) 



It follows immediately from these considerations that 
a boundary wave packet centered on these values must 
be correctly propagated by the vertex amplitude, if the 
vertex amplitude is to give the Einstein equations in the 
classical limit. 

The simplest wave packet we may consider is a diagonal 
gaussian wave packet 



where 



and 



(5) 



(6) 



(7) 



The normalization factors \/d~i = \J1i + 1, here and be- 
low, are required by the conventions we use in this paper, 
where the 3j-symbols are normalized. For details on the 
conventions, see [3]. The constants tr and Q are fixed by 
the requirement that the state is peaked on the value 
in — *o also when we change the pairing at the vertex. 
That is, by the requirement that all angles of the tetra- 
hedron are equally peaked on i„ ~ iq. It was shown in [7 
that this requirement fixes these constants to the values 
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(8) 



In other words, the state considered is formed by a gaus- 
sian state on the spins, with 0-phases given by the ex- 
trinsic curvature and by a "coherent tetrahedron" state 
(see H) 



(9) 



for each tetrahedron. 

Let us write the wave packet ([5]) as an "initial state" 
times a "final state" by viewing the process represented 



by the spacetime region described by the Regge cell as a 
process evolving four tetrahedra into one. That is, let us 
write this state in the form 



(10) 



where i'n — (ii, i4). Then we can test the classical limit 
of the vertex amplitude by computing the evolution of 
the four "incoming" tetrahedra generated by the vertex 
amplitude 

0(«) = X! ^0"™'4,*)V'»™t(jnm,4) (11) 



where i is 15, and comparing with '(/'(*)• If the vertex 
amplitude has general relativity as its classical limit, then 
we expect that in the semiclassical limit, namely for large 
Jo, the evolution should evolve the "initial" boundary 
state 'ipiijnrm'i'n) i^to a final state 0(i) which is still a 
wave packet centered on the same classical tetrahedron 
as the state ^p{i) given in ([9|. That is, 0(i) must be a 
state "similar" to plus perhaps quantum corrections 
representing the quantum spread of the wave packet. 

We have tested this hypothesis numerically, under a 
drastic approximation: replacing the gaussian depen- 
dence on the spins with a state concentrated on jnm = jo- 
That is, we have tested the hypothesis in the t 00 
limit. We are not sure this approximation makes sense, 
because in this limit the G phase dependence drops, and 
this may jeopardize the semiclassical limit. This approxi- 
mation is therefore only a first tentative hypothesis, that 
allows us to perform numerical calculations. Numerical 
calculation have proven useful in quantum gravity [9J, 
and the numerical investigation of the new vertex has al- 
ready begun |10j : here we consider a different angle from 
which the problem can be explored. 

The hypothesis we want to test is thus the following. 
We want to compare the evolved state 

4 

(P{i)= J2 W{iu.:,l4,i)'[[Hln) (12) 
ii...i4 n—1 

with the coherent tetrahedron state Q, where 

W{in) :-M^(jnm,«n)|,„„=jo- (13) 

If the function (p{i) turns out to be sufficiently close to 
the coherent tetrahedron state '0(*), we can say that, un- 
der the hypotheses given, the flipped vertex amplitude 
appears to evolve four coherent tetrahedra into one co- 
herent tetrahedron, consistently with the flat solution of 
the classical Einstein equations. 

From |2j the flipped vertex reads in the present case 

where 

4 / . ■ \ 
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The 15j-symbol is the contraction of five 4-valent inter- 
twiners, each constructed with two Wigner 3j-symbols 
and each normahzed with a \fdi factor. The matrix in- 
dicates the Wigner 3j-synibol. The indices a„ are in the 
representation jo ^^nd the indices 5„ and c„ in the repre- 
sentation jo /2. The quantities iai...a4 are the components 
of the invariant tensors that define the intertwiner i. For 
technical details, see 

We have compared the two functions (coherent 
tetrahedron) and (^(i) (evolved state) for the cases jg = 2 
and jo = 4. The numerical results are shown in the fig- 
ures below. The overall relative amplitude of ^(«) and 
^(i) is freely adjusted by fixing the normalization con- 
stant N and therefore is not significant. The quantity 
*mcan is the mean value of It gives the position of 
the wave packet. The quantity cr/2 is the corresponding 
variance. It gives the (half) width of the wave packet. 
In Fig.l and Fig. 3 we compare the modulus square of 
the wave function (for the two values of jo). In Fig. 2 
and Fig.4 we compare the modulus square of the discrete 
Fourier transform of the wave function: n stands for the 
n"' multiple of the fundamental frequency 27r/jo. 
Case jo = 2 : 
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FIG. 1: jo = 2. Modulus square of the amplitude. Left: co- 
herent tetrahedron (Zmeaii±f7"/2 = 2.54±0.39j. Right: Evolved 
state (iraean±(j/2 = 2.54±0.46;. 




FIG. 2: jo = 2. Modulus square of the (discrete) Fourier 
transform of the amplitude. Left: coherent tetrahe- 
dron (nmean^(y/2 — 1.25±0.27j. Right: Evolved state 
fn„,a„±o-/2 = 1.15±0.3i;. 



The agreement between the evolved state and the co- 
herent tetrahedron state is quite good. Besides the over- 
all shape of the state, notice the concordance of the mean 
values and the widths of the wave packet. Considering 
the small value of jo, which is far from the large scale 
limit, and the r — s- oo limit we have taken, we find this 
quite surprising. 

The same pattern repeats in the jo = 4 case: 
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FIG. 3: jo = 4. Modulus square of the amplitude. Left: co- 
herent tetrahedron (imean^o- /2 = 4.88±0.56/ Right: Evolved 
state (i,nean±cr/2 = 4.85±0.96;. 
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FIG. 4- jo = 4. Modulus square of the (discrete) Fourier trans- 
form of the amplitude. Left: coherent tetrahedron (n,neani 
a/2 = 2.25 ±0.32/ Right: Evolved state fn„<.an±o-/2 = 
2.08±0.59j. 

The technique developed here is complementary to the 
one developed in [11 , which is based on computing n- 
point functions. We expect that this technique could be 
better developed. 

The numerical results presented above are preliminary 
and tentative, but they appear to support the expecta- 
tion that the flipped vertex might in fact give general 
relativity in the classical limit. 

We thank Eugenio Bianchi and Etera Livine for numer- 
ous useful exchanges and suggestions. 
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